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ABSTRACT 

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferom- 
eter under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. 
We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50° diameter 
fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the 
MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements 
of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, 
and used it to produce ~ 15' angular resolution maps of the two fields. We perform a blind source 
extraction using these confusion-limited images, and detect 655 sources at high significance with an 
additional 871 lower significance source candidates. We compare these sources with existing low- 
frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis 
algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency 
surveys in these regions of the sky and with statistical distributions of point sources derived from 
Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 
110 MHz to 200 MHz band. 

Subject headings: Surveys - dark ages, reionization, first stars - Instrumentation: interferometers - 
Methods: data analysis 



l introduction 

The study of the origin and evolution of the Universe 
draws upon observations of phenomena at a large range 
of distances and look-back times, connecting the ini- 
tial conditions probed by the cosmic microwave back- 
ground to present-day conditions dominated by galaxies 
and clusters of galaxies. A major chapter of this his- 
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tory has yet to be examined — the chapter that corre- 
sponds to redshifts from z = 1000 to z = 6 and comprises 
the Dark Ages and the Epoch of Reionization (EoR) of 
the Universe. The EoR, in particular, marks a major 
milestone when the first stars and galaxies formed and 
reionized the intergalactic medium. It was recognized 
some time ago that studies of the redshiftcd 21-cm ra- 
dio emission from ne utral hydrogen would be a promis- 
ing pr obe of the EoR (iHogan fc Reesl ll979l: iScott fc Reei 
119901: iMadau et all 119971) . Indeed, the possible exis- 
tence of extensive regions containing significantly large 
amounts of neutral hydrogen at redshifts of z = 15 to 
z = 8 motivates an interest in developing highly sen- 
sitive low-frequency radio telescopes in order to detect 
the redshifted 21-cm signal. Not long ago it was recog- 
nized that a statistical detection of the patchy neutral 
hydrogen distribution during the EoR should be pos- 
sible with an array w ith a modest c ollecting area and 
a large field of view (jMorales fc Hew itt 2004). In the 
past decade considerable efforts have been made towards 
this goal, through advances in theoretical modeling of 
the EoR signature (see e.g. iPritchard fc Loebl 1201 IL 
[Morales fc Wvithd [20Tol and iFurlanetto et al.ll2006r for 
recent reviews), as well as through the development 
of new instrumental approaches to measure the red- 
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shifted 21-cni sign a l (se e e.g. i Bowman fc Roger s! 120101 
Chippendali 120091 iLqi isdalc ct all 1200^. Tinga v et al] 



in prep.l Parsons et al.ll2QlQnRottgering et al.ii2Q06l and 



Paciga et a l. 2011) 



One particularly daunting challenge for these experi- 
ments is emission from foreground astrophysical sources, 
which is at least two to three orders of magnitude 
brighter than the redshifted 21-cm signal (in the total 
intensity as well as in the magnitude of the spatial fluctu- 
ations, see e.g. Shaver et al. 1999: 'Bcrnardi et al.l 120091 : 
Een et al. 2009; dc Olivcira-Costa ct al. 2008). The fore- 
ground emission arises mainly from synchrotron and free- 
free processes, and therefore is highly likely to have 
an intrinsically smooth radio spectrum. The 21-cm 
signal, however, is likely to be produced under condi- 
tions that vary rapidly with redshift, and, if this is the 
case, will appear to have rather sharp spectral features. 
Techniques have been developed to exploit this differing 
spectral behavior in order t o separate and remove the 
foreground contamination (iFurlanetto fc Briggsl 120041: 
McQuinn ct al. 2006: Bow man et al.ll2009l: iHarker et ail 
2009t iLiu et al.ll2009t iLiu fc Tegmarkll201ir ). However. 



any calibration imperfection or instrumental defect has 
the potential to introduce distortions into measured spec- 
tra, and thereby to mix the extremely bright foreground 
emission with the signal from the redshifte d 21-cm sig- 
nal in ways that are difhcult to disentangle (|Datta et al.l 
120101 ). Thus, an equally daunting challenge is to learn 
how to calibrate any new instrumentation that is being 
developed for these observations with extremely high fi- 
delity. 

The Murchison Widefield Array (MWA, 
iLonsdale et al.l 120091 : iTingav et al.l lin prep.f ) is a 
new array being constructed to characterize the 21-cm 
signal during the EoRFI. In addition to the study of 
the EoR, other key science goals of the MWA include 
the study of radio transients, the study of the helio- 
sphere and ionosphere, and low-frequency Galactic and 
extragalactic studies. These four goals and potentially 
others are addressed by an array made of a large number 
of small antenna elements that simultaneously give a 
large collecting area and a large field of view. This 
is a departure in many ways from a traditional radio 
array design, with phased arrays of dipoles constituting 
the fundamental antenna elements, digitization early in 
the data stream, and full correlation of a large number 
of baselines. This design promises large improvements 
with regard to wide-field surveys and detection of the 
EoR, but, at the same time, it poses new challenges, 
especially with regard to calibration and imaging of 
the large field of view and compensation for the effects 
of the ionosphere. The instrument is currently under 
construction; work on a 128-element array commenced 
in early 2012 at the radio-quiet Murchison Radio- 
astronomy Observatory (MRO) in Western Australia. 
As a first step in demonstrating the new technologies 
required for the MWA, a 32-element prototype was 
built at the MRO site prior to the build-out to 128 
elements. This prototype system was operated for 
two years in campaign mode, and underwent a cycle 
of equipment installation, testing, and redesign as 

^"^ Additionally, the MWA has been designated an official pre- 
cursor instrument for the Square Kilometre Array (SKA). 



necessary. Beginning in March 2010, the prototype was 
used for initial science o bservations, a nd ha s already 
yielded several results. lOberoi et al.l ()2011[ ) presents 
findi ngs from an inve stigation of solar radio emission, 
and Ord e t aD (j2010( ) present wide-field images using 
a prototype real-time imaging and calibration pipeline. 
Herein, we also report results based on data obtained 
during this initial science run. 

The goals of the measurements and analysis presented 
here are to verify the performance of the MWA subsys- 
tems and the 32-element prototype array, to explore tech- 
niques for future EoR experiments, and to deepen our 
understanding of the astronomical foregrounds. We ob- 
served two overlapping fields at high Galactic latitude, 
each 50° across. One field was identified for possible 
EoR studies in the future, and the other was chosen to 
have a very bright radio source at its center to facili- 
tate calibration. The observations are deep in the sense 
that they combine data from a large range of hour angles 
and multiple snapshot images to improve sensitivity and 
image fidelity; developing such techniques for deep ob- 
servations is critical for the success of future EoR experi- 
ments and other scientific investigations with the MWA. 
We have developed a data reduction pipeline that im- 
plements wide-field calibration and imaging algorithms, 
compensates for the direction-dependent and changing 
primary beam as different snapshots are combined, and 
automatically extracts sources from the images. We com- 
pare our results directly to the results from other sensi- 
tive low-frequency radio surveys in the southern hemi- 
sphere, and we compare them statistically to the results 
from surveys carried out in the northern hemisphere. We 
use these comparisons to assess the performance of the 
MWA prototype and the wide-field imaging and calibra- 
tion algorithms. We assess the completeness and reli- 
ability of our point source catalog through comparison 
to these surveys. We make a number of simplifying as- 
sumptions in this first phase of analysis; future work will 
refine the techniques and algorithms until the stringent 
calibration requirements of EoR experiments with the 
full MWA can be met. 

2. LOW-FREQUENCY RADIO SURVEYS 

In this work we make extensive use of the results of 
previous sensitive low-frequency radio surveys to verify 
the performance of the MWA and to provide external 
data for calibration. We summarize here the properties 
of the surveys used in our comparisons. 

Th e Molonglo Reference Catalog (MRC; iLarge et al.l 
I1981D is the product of a blind survey at 408 MHz 
that covered nearly all of the southern sky to moder- 
ate depth. The catalog covers all right ascensions in 
the declination range from —85.0° to -1-18.5°, excluding 
the area within 3° of the Galactic plane. The obser- 
vations were conducted with the Molonglo Radio Tele- 
scope in a 2.5 MHz wide band with a synthesized beam 
of 2.62' X 2.86'sec(<5 + 35.5°) in width. The MRC has 
a stated completeness limit of 1 Jy, although it contains 
sources down to a flux of ^0.7 Jy. 

At frequencies below 408 MHz, there have been many 
targeted observations of known sou rces in the southern 
sky. The Culgoora Circular Array ()Sleelll995l) was used 
to observe Galactic and extragalactic sources selected 
from existing higher frequency surveys. The observa- 



MWA-32T Low Frequency Imaging 



tions were made at frequencies of 80 MHz and 160 MHz. 
The beam size was 3.70' x 3.70'sec(5 + 30.3°) at 80 MHz, 
and 1.85' x 1.85' sec(5 + 30.3°) at 160 MHz. The Umiting 
flux density was 4 Jy at 80 MHz, and 2 Jy at 160 MHz. 
However, only a small patch of sky around each selected 
source was imaged. Although flux densit y uncertain - 
ties are not directly reported in this list, iSled ()1977f ) 
note that the standard deviation in the flux density for 
sources measured with the Culgoora array is ~ 13% 
for the brightest sources, and ~ 39% for the faintest 
sources, with a pote ntial systemat i c flux scale depression 
of ^ 10%. Similarlv. iJacobs et al.l (|2011l ) present results 
from PAPER, an array of east-west polarized dipoles, 
that were obtained over the 110 MHz to 180 MHz band. 
The results were derived from multi-frequency synthesis 
maps of the entire sky south of a declination of 10° hav- 
ing a resolution of 26'. A sample of 480 sources with 
fluxes greater than 4 Jy in the MR C were identif e d and 
measured in these PAPER maps. iJacobs et al.l ()2011[ ) 
find a 50% standard deviation in their fluxes relative to 
values obtained from the MRC and Culgoora source lists. 
They quote a flux limit of 10 Jy for the sources in their 
catalog. 

There are several on going efforts t o perform low- 
frequency blind surveys. iPandevI (|2OO60 presents results 
from a survey that used the Mauritius Radio Telescope 
(MRT) to image ~1 steradian of the sky at 151 MHz 
and to thereby produce a catalog of 2782 sourceJ*^. 
The deconvolved images achieve an angular resolu- 
tion of 4' X 4.6'sec(5 -I- 20?14) and a root-mean-square 
(RMS) noise level of approximately 300 mJy beam" 
^avak et al.' "M^. The TIER GMRT Sky Survey 
(TGSS, Sirothia et all lMTI ) is producing a 150 MHz sur- 
vey of the sky at declinations above —30°. Each pointing 
covers ~ 7 square degrees and yields a map that reaches 
an RMS noise of ^ 8 mJy beam" at an angular reso- 
lution of ~ 20". The flux density scales of the maps are 
limited by systematic errors and have relative errors of 
25%. As of 2012 January, the TGSS website^ reports re- 
sults from images of approximately 2600 square degrees 
of the southern sky. 

Surveys which primarily cover the northern sky have 
also been carried out at low frequencies. The most ex- 
tensive wide-field uniform survey near our observi ng fre- 
quen c y of 150 MHz is the 6C survey (Baldwin et al.l 
[19851: iHaleTeFall [1988L ITQQOI HqML ll993blla] V The 6C 
survey covered the northern sky above declination 30° 
with a sensitivity of 200 niJy; the angular resolutio n was 
4.2' X 4.2' CSC (5. The 7C survey (|Hales et al.ll2007t ) cov- 
ers 1.7 sr of the northern sky to a greater depth and at 
higher resolution than the 6C survey. We have chosen 
to use the 6C survey for comparison to our results in 
this paper because it covers a somewhat larger sky area, 
and has served as the basis for other investigations of 
EoR foregrounds (in particular, . Pi Mattco et al. 2003)- 
At lower frequencies, iCohen etal.1 (|2007D have used the 
VLA 74 MHz system to perform a survey of the sky north 
of declination S = —30°. This survey, known as the VLA 
Low- Frequency Sky Survey (VLSS), produces maps with 
an 80" angular resolution which achieve a typical RMS 



noise level of 100 mJy beam" . Each VLSS image is 
14° X 14° across in order to fully image the VLA primary 
beam, which has a FWHM diameter of 11.9° degrees. 
ICohen et al.l (|2007[ ) perform a blind source extraction on 
the VLSS maps and produce a source catalog of 68,311 
radio sources above a significance level of 5cr. They quote 
a 50% point source detection limit of 0.7 Jy beam" . 

For completeness we note that additional low- 
frequency surveys include the Sydney University Mo- 
, . „, „ ,„.„ ,,„ " ' ^• |T99l . the 

19971) . and 



longlo Sky Survey (843 MHz, Bo ck et al 

Miyun survey (232 MHz, Zhang et al.l 

the Levedev Physical Inst itute Survey (102.5 MHz, 

iDagkesamanskii et al.l I2000D . We have not extensively 

compared our results to results from these surveys. 

Low frequency surveys have also been used to 
study t he distribution of spe ctral indices of radio 
sources. iDe Breuck et al.l (|2000) used results from the 
MRC and f rom the Parkes-MI T-NRAO 4.85 GHz sur- 
vev (PMN. IWright et al.lll99llG rifiith fc Wrightlll993l: 
Griffith et al.l 1199^ 'Condon et al.i il993t iTasker et al.l 



° Electronic c atalogs are available 

|http : / /www . rri . res . in/surveys /MRT| 

^" |http : //tgss . ncra . t if r . res ■ in/150MHz/tgss .htmTj 
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1994; Grif fith et al.i 1995; Wrig ht et al.'ril996( ) to st"udv 
the distribution of spectral indices of sources in the 
southern sky. They also carried out similar comparisons 
of the results from the Westerbork Northern Sky Sur- 
vey (325 MHz; 'Rengelink et al."1997[) and the Texas Sur- 
vey (365 MHz; Douglas et al. 199^) with results from the 
NRAO VLA Sky Survey (1.4 GHz; ICondon et allfl^ 
in the northern hemisphere. The spectral index distri- 
butions showed significant differences between samples 
selected at low frequencies and samples selected at high 
frequencies. The combined MRC-PMN source list was 
also used to generate a sample of ultra-steep spectrum 
sources. 

In Section 17.11 we carry out source-by-source compar- 
isons of our survey results to those of the MRC, the Cul- 
goora fiux density measurements, the PAPER flux den- 
sity measurements, and the TGSS. There is no overlap 
at present between our survey and the Mauritius survey, 
but comparisons should become possible when the anal- 
ysis of the Mauritius data is completed. We also carry 
out statistical comparisons of our survey results to those 
of the 6C survey by com paring source counts, and to the 
IDe Breuck et al.l ()2000D spectral index catalogs by com- 
paring spectral index distributions. 

3. THE MWA-32T INSTRUMENT 

The Murchison Widefield Array 32-Tile prototype 
(MWA-32T) was built and operated for the purpose of 
verifying the performance of MWA subsystems in prepa- 
ration for building a larger, more capable array. As noted 
above, construction of a 128-tile array has commenced 
and is expected to be complete later this year (2012). 
We summarize t he de sign here; the reade r is referr ed to 
iLonsdale et al.l ()2009[ ) and iTingav et al.l (|in prep. I ) for 
more detailed descriptions. 

The MWA-32T was designed to cover a frequency 
range from 80 MHz to 300 MHz, with an instanta- 
neous bandwidth of 30.72 MHz at a spectral resolution 
of 40 kHz. The array consisted of 32 antenna tiles, 
which served as the primary collecting elements of the 
array. The tiles are designed to have an effective col- 
lecting area larger than 10 m^ in the MWA frequency 
band, and to provide a steerable beam which can be 
pointed up to 60° from zenith while maintaining a sys- 
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tern temperature which is dominated by sky noise within 
the MWA band. The primary beam of the tiles is fre- 
quency and position dependent, with a FWHM size at 
zenith of roughly 25°/{iy/150 MHz). Each tile consists 
of 16 dual-polarization, active dipole antennas laid out 
over a metal mesh ground screen in a 4 x 4 grid with a 
1.1 meter center-to-center spacing. Each dipole antenna 
consists of vertical bowtie elements that feed a pair of 
integrated low-noise amplifiers (LNAs) located within a 
tube at the juncture of the orthogonal arms of the dipole. 
The antennas are designed to have low horizon gain to 
reduce terrestrial RFI contamination, and to have a low 
manufacturing cost. 

The signals for the two polarizations are processed in 
parallel. For each polarization, the signals from the 16 
dipoles on each tile are carried over coaxial cable to an 
analog beamformer, where they are coherently summed 
to form a single tile beam. A system of switchable ana- 
log delay lines is used to apply an independent time delay 
to each of the dipole signals, allowing the tile beam to 
be steered on the sky. The delay lines employ a series 
of 5 switchable traces, each differing by a factor of two 
in length, with the shortest trace introducing a nominal 
delay of 435 ps. This allows for 32 discrete delay set- 
tings for each of the input signals. The discretization of 
the delays implies that the primary beam can only be 
steered in discrete steps, and so can only coarsely track 
a sky field. The summed signal is amplified and sent over 
coaxial cables to the MWA digital receiver for digitiza- 
tion. 

Each MWA digital receiver node services 8 tiles. The 
16 received signals are first subjected to additional fil- 
tering and signal conditioning for low-frequency rejec- 
tion, anti-aliasing, and level adjustment. The signals are 
then digitized at baseband by eight dual 8-bit analog-to- 
digital converter (ADC) chips operating at a sampling 
rate of 655.36 MHz. The data stream from each ADC is 
fed to a digital polyphase filterbank (PFB) implemented 
in FPGA hardware which produces 256 frequency chan- 
nels, each 1.28 MHz wide. A sub-selection of 24 of these 
channels (a bandwidth of 30.72 MHz) are transmitted 
via optical fiber to the correlator. 

At the correlator, the data streams from each receiver 
are processed by a second stage PFB to obtain a fre- 
quency resolution of 10 kHz. The signals are then cross- 
multiplied to produce a 3072 channel complex spectrum 
for each of the 2080 correlation products. These com- 
prise the four polarization products for all pairs of tiles 
as well as the autocorrelations. The visibilities are aver- 
aged into 40 kHz wide channels and integrated for 50 ms 
due to output data rate constraints. During the 32T ob- 
serving campaign described in this work, the correlator 
was operating at a 50% time duty cycle due to hardware 
limitations. The visibilities are captured and averaged 
into 1 second integrations before being written to disk. 

4. OBSERVATIONS 

Observations were conducted with the MWA-32T in 
March 2010 during a two-week campaign (X13) when 
personnel were present on-site to operate the instrument. 
Data were taken in three 30.72 MHz sub-bands centered 
at 123.52 MHz, 154.24 MHz, and 184.96 MHz in order 
to give (nearly) continuous frequency coverage between 
~110 MHz and ~200 MHz. During the observations. 



TABLE 1 
Journal of Observations 



Field 


Frequency 


Date 


Number 


Observation'' 




(MHz) 




of Scans 


Time (minutes) 


EoR2 


123.52 


2010/03/24 


35 


208 






2010/03/28 


26 


154 




154.24 


2010/03/22 


35 


208 






2010/03/26 


18 


107 






2010/03/28 


10 


59 




184.96 


2010/03/21 


35 


208 






2010/03/25 


35 


208 






2010/03/26 


18 


106 






2010/03/29 


36 


214 


Hydra A 


123.52 


2010/03/24 


8 


39 






2010/03/28 


7 


34 




154.24 


2010/03/22 


8 


39 






2010/03/26 


5 


24 






2010/03/28 


3 


15 




184.96 


2010/03/21 


8 


39 






2010/03/25 


8 


39 






2010/03/26 


5 


24 






2010/03/29 


9 


44 



'^ The effective integration time is less than half of this observation 
time due to the 50% duty cycle of the correlator and additional 
flagging. 

the beamformers were used to steer the beam in steps 
as the fields crossed the sky. This stepped steering is 
a consequence of the discretization of the analog delay 
lines in the beamformer. The typical sequence was to 
steer the beam to a new position, observe at a particular 
frequency for five minutes (without tracking), and then 
steer the beam again. Thus, the measurements can be 
considered to be a series of short drift scans. 

The observing time was divided between two fields. 
One field was centered on the bright extragalactic source 
Hydra A at RA(J2000) = 9^ 18"" 6", Dec.(J2000) = 
— 12° 5' 45" to facilitate calibration. The other covered 
the EoR2 field, centered at RA(J2000) = 10^ 20" 0", 
Dec.(J2000) = -10° 0' 0". The EoR2 field is one of 
two fields at high Galactic latitude that have been iden- 
tified by the MWA collaboration as targets for future 
EoR experiments. It also had the advantage of being 
above the horizon at night during the observing cam- 
paign. Although the centers of the Hydra A and EoR2 
fields are separated by 15.3°, there is considerable over- 
lap between them since the half power beam width of 
the primary beam is ^ 25° at 150 MHz. A total of 61 
~5 minute scans of the Hydra A field and 248 scans of 
the EoR2 field were obtained in interleaved sequences 
over the course of the observing sessions. Table [1] gives 
a journal of the observations. 

5. DATA REDUCTION STRATEGY 

5.1. Instrumental Gain Calibration 

The MWA antenna tile architecture poses several non- 
traditional calibration issues due to both the nature of 
the primary beam and the wide field-of-view. The pri- 
mary beam is formed by the summation of beamformer- 
delayed zenith-centered dipole responses. The beam- 
former delays are periodically changed to track a field 
across the sky. Although this moves the center of the 
primary beam as intended, the overall shape of the beam 
changes as well. For a given set of beamformer delays, 
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the beam is fixed relative to the tile, and therefore moves 
relative to the sky as the Earth rotates. As a further com- 
plication, the time and direction dependence of the pri- 
mary beam response is different for the two polarizations 
of the crossed dipoles in the array. This leads to appar- 
ent polarization in inherently unpolarized sources due to 
the different responses in the two orthogonal dipole po- 
larizations unless the appropriate corrections are applied 
in the analysis procedure. 

Methods for measuring and calibrating the primary 
beam have been dev eloped for the full MWA system by 
iMitchell et al.l (|2008| ) , who plan to use a real-time system 
(RTS) to calibrate and image MWA data. Their method 
performs a calibration of the instrument in real time, 
solving for direction- and frequency-dependent gains for 
each antenna based on the simultaneous measurement of 
multiple known bright sources across the field of view. 
This method was developed for use in a 512-tile array, 
where the instantaneous sensitivity and ww-plane cover- 
age enable the measure ment of several hundre d sources in 
each 8 -second iteration ()Mitchell et al.ll2008D . lOrd et al.l 
(|2010( ) have successfully demonstrated the use of a mod- 
ified version of the RTS in order to calibrate and image 
data from the MWA-32T. However, the reduced sensi- 
tivity of the MWA-32T array makes this full calibration 
challenging. For the MWA-32T system, the data rate 
is sufficiently low that real-time calibration and imag- 
ing are not necessary, and the raw visibility data can be 
captured and stored. We have chosen to pursue an alter- 
native data reduction pipeline based on more traditional 
calibration and imaging software which allows us to use 
the full visibility data set in order to perform a detailed 
investigation of the calibration and imaging performance 
of the MWA-32T instrument. 

Without the ability to directly measure the primary 
beam for each tile, we instead assume a model and use 
it to account for the instrumental-gain direction depen- 
dences. Knowledge of the primary beams is also needed 
for an optimized weighting of the maps when combining 
them to obtain deeper maps (see Section [6^ . It is likely 
that precise characterization of individual tile beams will 
be necessary to achieve dynamic range sufficient for accu- 
rate foreground subtraction and EoR detection, but this 
is not attempted in the work described herein. 

For present purposes, we assume the polarized pri- 
mary beam patterns are identical across all tiles, and can 
be modeled by simply summing together the direction- 
dependent complex gains of the individual dipoles in the 
tile, i.e., mutual coupling between elements and tile-to- 
tile differences are ignored. We model the complex beam 
patterns of an isolated individual dipole for both the 
north-south (Y) and east-west (X) electric field polariza- 
tions using the WIPL-D Prcj^ electromagnetic modeling 
software package. A tile beam pattern is then computed 
by summing the 16 dipole responses with the dipoles as- 
sumed to be at their nominal locations in a tile and with 
the individual responses modified by the nominal ampli- 
tudes and phases introduced by the beamformer for the 
given delay settings. Figure [T] shows cuts through power 
patterns (square modulus of the complex beam) at zenith 
and at a pointing direction 28° east of zenith for both the 
X and Y polarizations. Model beam patterns were cal- 

^^ http://www.wipl-d.coirL 
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Fig. 1. — West-east cuts through simulated MWA antenna tile 
patterns for a zenith (a,b) and a 28° easterly pointing (c,d) at 
150 MHz. Panels (a) and (c) show the X polarization dipole power 
response pattern, while panels (b) and (d) show the Y polarization 
power response pattern. The scale is logarithmic, with arbitrarily 
normalized decibel units. The polarization-dependent gain struc- 
ture is clearly visible in the sidelobes. 



culated at frequency intervals of 2 MHz, since they vary 
significantly across the MWA frequency band. 

We assume that this model fully describes the direc- 
tion dependence of each tile. We do, however, allow for 
a different overall, i.e., direction- independent, complex 
gain for each til e. We follow the Jones matrix formalism 
as presented bv iHamaker et all ()1996[) . The instrumen- 
tal model then takes the following form for a single tile 
at a single frequency: 



where 



VA = GaBga, 



eA 



(1) 

(2) 



is the incident electric field at tile A, decomposed into 
linear E-W and N-S polarizations. 



VA 



is the vector of measured antenna voltages, 



G. 



5A,x 







9A.y 



(3) 



(4) 



is the matrix of direction-independent complex gains for 
an antenna, and 



B = 







(5) 



is the matrix of direction- and frequency-dependent but 
tile-independent gains due to the primary beam shape 
(we represent spatial coordinates with 9 and <l>). We ne- 
glect the feed-error "D" matrix of lHamaker et al] ([1996); 
in other words we assume that the sensitivity of the X- 
polarization response of a dipole to Y-polarization radia- 
tion (the cross-polarization) is zero, and vice versa. This 
is likely to be a good approximation since ideal dipoles 
have zero cross-polarization by definition. In reality, var- 
ious effects, such as the finite thickness of our dipole 
elements, interactions between structures in neighbor- 
ing dipoles, or projection effects may produce a nonzero 
cross-polarization response. In this paper we restrict our 
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imaging and analysis to these two senses of linear polar- 
ization and their combination as total intensity. Errors 
caused by neglecting cross-polarization effects are second 
order in the small off-diagonal elements of the D-matrix. 

The strategy we adopt for the data reduction is first to 
analyze short snapshots wherein the settings in the ana- 
log beamformer were static so that the primary beam 
pattern can be taken as constant over the duration of 
each snapshot, and any gain changes due to the side- 
real motion of the sky relative to the beam can be ne- 
glected. In this regime, the direction-dependent gain can 
be factored out of the response and corrected in the image 
plane in the resulting map. Under this approximation, 
we are able to use standard tools for radio astronom- 
ical data analysis for much of the processing. Finally, 
the frequency dependences of the overall antenna com- 
plex gains are in principle determined by modeling the 
summed spectra of the bright sources in the field used in 
the calibration. 

Standard calibration procedures rely on being able to 
observe a field containing a strong source with easily 
modeled structure that substantially dominates the vis- 
ibilities. For the two fields presented in this work, Hy- 
dra A is the strongest source in the field. iLane et aD 
(|2004f) present low-frequency images that show that while 
it is quite extended at the VLA's resolution, most of the 
flux is contained within a region that is a few arcminutes 
in radius. Since this extent is smaller than the angular 
resolution of the MWA-32T, we were able to treat Hy- 
dra A as a point source in our calibration analysis. One 
might expect that for a large field-of-view instrument 
such as the MWA, we would also need to include several 
or even many additional strong sources in the calibra- 
tion model with known direction-dependent gains. We 
therefore experimented with calibration models that in- 
cluded several point sources in addition to Hydra A, but 
we found that the complex tile gain solutions were not 
significantly changed. We therefore simply used Hydra A 
as the only calibration source in subsequent analysis. It 
should be noted that this is a potential source of error. 

5.2. Ionosphere 

At the low radio frequencies of the MWA, position- and 
time-dependent variations in the electron density of the 
ionosphere cause variations in propagation times, which 
appear in the visibility data as frequency- and time- 
dependent phase shifts. For the short baselines of the 
MWA-32T, these variations are, except at times of ex- 
treme ionosphere disturbance, refractive in effect, i.e., 
they simply cause apparent changes in the positions of 
point sources on the sky. These positions shifts may be 
different in different directions, especially over a wide 
field such as that of the MWA, and consequently may 
lead to distortions in the derived images. 

Ionospheric effects have been quantif ied by studies 
with o ther low-frequency interferometers. [Baldwin et al.l 
()1985[) . using the Cambridge Low- Frequency Synthesis 
Telescope at 150 MHz, found that the ionosphere typi- 
cally caused 5° (RMS) phase variations on 1 km baselines 
on sub-day timescales, which were uncorrelated from day 
to day. They also remarked that ionospheric irregular- 
ities on large spatial scales, most likely related to the 
day-night cycle and strongly correlated from one day to 
the next, could induce apparent position shifts of sources 



of up to 20'. iKassim et all ()2007[ ) found, with the VLA 
operating at 74 MHz, that during times of moderate iono- 
spheric disturbance relative position variations on short 
times cales across a 2 5 ° field of view were at most 2'. Sim- 
ilarly, iParsons et al.l (|2010t ) , using observations of bright 
sources with the PAPER array at 150 MHz, found short- 
term small (typically less than 1') position offsets that 
were not correlated from day to day, and long-term large 
(up to 15') position offsets which were correlated from 
day to day, and were mainly in the zenith direction. 

In our analysis, we average snapshots taken with the 
center of the field within several hours of the local merid- 
ian over a period of eight days. The results in the papers 
cited above suggest that uncorrelated short-term vari- 
ations in source directions will be significantly smaller 
than our beam size of ~ 15' and, furthermore, that they 
will tend to average out when images derived from in- 
dividual snapshots are combined. We therefore neglect 
them. Long-term correlated variations in source direc- 
tions may be comparable to our beam size, and may not 
average to zero as we combine snapshot images. How- 
ever, our calibration strategy, described in Section 15.11 
will tend to remove any ionospheric offset at the position 
of Hydra A through the phase terms in the direction- 
independent gain solutions. It is possible, depending on 
the behavior of the ionosphere during the present ob- 
servations, that long-term differential position shifts of 
several arc minutes might be present in our final images. 
We investigate this possibility through comparisons of 
the positions of our extracted sources with source posi- 
tions listed in published catalogues. 

Our neglect of short-term ionospheric effects is justi- 
fied only because of the small baselines (<350 m) of the 
MWA-32T array. For the longer baselines (~3 km) of 
the full MWA, we believe that these effects will need to 
be corrected to achieve the dynamic range required for 
many of the science goals. 

6. REDUCTION PIPELINE 
6.1. Initial Processing and Editing 

We developed a calibration and imaging pipeline based 
on the NRAO Common Astronomical Software Appli- 
cationsPI (CASA) package and additional tools that we 
developed in Python and IDL. The pipeline uses a se- 
ries of short observations to generate "snapshot" images 
which are weighted, combined, and jointly deconvolved 
to produce final integrated maps. 

In the first stage of the pipeline, the visibilities were 
averaged over 4-s intervals and converted from the MWA 
instrumental format into UVFITS files. The MWA-32T 
correlator does not perform fringe-stopping (the correla- 
tion phase center is always at zenith) , so phase rotations 
were applied to the visibilities to track the desired phase 
center. As a part of this process, data corrupted by RFI 
were flagged for later exclusion from the analysis. The 
data were then imported into CASA. Additional editing 
was done to flag data affected by known instrumental 
problems or RFI. Approximately 25% of the data were 
flagged at this stage, mainly because a problem in the 
data capture software corrupted 480 of the 2080 correla- 
tion products. 

^^ http://casa.nrao.edu 
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6.2. Calibration of Antenna Gains 

Calibration was performed separately for each snap- 
shot with CASA, using Hydra A as the gain and phase 
reference (as discussed in Section 15. ip . Although Hy- 
dra A was not at the center of the primary beam during 
observations of the EoR2 field, it was still strong enough 
to substantially dominate the visibilities. Model visibil- 
ities were calculated using a point source model for Hy- 
dra A, assuming an unpolarized flux of unity. The over- 
all frequency-dependent flux scale of the data was set at 
a later point, along with a correction for the direction- 
dependent gain. Time- independent channel- by- channel 
gain factors were calculated for each tile using the task 
bandpass. After this was done and the visibilities cor- 
rected for the gain factors, time-dependent overall tile 
gain factors were calculated on a 32 second cadence us- 
ing the task gaincal. The factors determined in the two 
tasks give the frequency- and time-dependent ^a.x and 
5A,y terms of Equation |4l 

The (7a, X and gA.y terms were examined for temporal 
stability and spectral smoothness; regions where devia- 
tions were apparent were flagged. Such deviations were 
rare, and an important outcome of this analysis is the 
recognition that the MWA antenna gains are quite stable 
over frequency and time. In fact, the gains even tended 
to be stable from one day to the next. However, some 
complicated variations in gain as a function of frequency 
were identified. These were associated with damaged ca- 
bles and connectors that have since been replaced or re- 
paired. 

6.3. Snapshot Imaging 

The data from each snapshot were subdivided into 
7.68 MHz wide frequency bands, and multifrequcncy syn- 
thesis imaging was performed for each snapshot using the 
CASA task clean. Images were made with a 3' cell size 
over a ~ 51°x ^ 51° patch of sky in order to cover 
the majority of the main lobe of the primary beam. 
The "XX" and "YY" polarizations were imaged sepa- 
rately. Conversion to the standard Stokes parameters 
was not performed at this stage, since, as discussed in 
Section 15.11 the gains for the two polarizations have dif- 
ferent direction dep endence. The "w-projection" algo- 
rithm (jCornwell et al. 2008) was used to correct for wide 
field-of-view effects, and to produce an image with an 
approximately invariant point spread function in each of 
the snapshot images. The images were dec onvolved usin; 
the Cotton-Schwab CLEAN algorithm feee lSchw"ablll98- 
down to a threshold of 1% of the peak flux in the image. 

A position-dependent "noise" map was also computed 
for each of the 7.68 MHz wide snapshot images by select- 
ing a 64-pixel by 64-pixel window around each pixel in 
the image and fitting a Gaussian to the central 80% of 
a histogram of the pixel values. This procedure was em- 
ployed because of the high point-source density in these 
maps. Throughout much of the area of these maps, it was 
impossible to identify a source free region from which 
to estimate the background noise fluctuations, and the 
presence of sources artificially skewed the noise estimates 
calculated strictly as the root-mean-square (RMS) of the 
pixel values. We found that this clipped histogram fit- 
ting procedure provided a more robust estimate of the 
RMS of the background noise. These noise maps were 



smoothed on a 1° scale to remove local anomalies intro- 
duced by extended or clustered sources. However, despite 
these procedures, some areas, particularly in especially 
crowded regions, still had anomalously high noise esti- 
mates. 

As discussed above, each snapshot was only ^5 min- 
utes in duration, and was obtained while the delay line 
settings in the analog beamformers were fixed. This al- 
lowed us to model the primary beam pattern of each 
tile as fixed relative to the sky for the duration of the 
snapshot. Our calculated model beams for each po- 
larization formed the frequency-dependent bx{0^(j)) and 
by (9,(1)) terms of Equation [51 These terms are time- 
dependent only in that they are different for each snap- 
shot. 

6.4. Snapshot Combination and Joint Deconvolution 

Deeper images were obtained by combining snapshot 
maps from a particular 7.68 MHz wide band according 
to: 
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(6) 



where /dirty is the integrated, primary-beam corrected, 
dirty map, the snapshots are distinguished by index z, 
Di is a snapshot dirty map, Bi is the primary beam cal- 
culated for that snapshot, and Ui is the fitted rms noise 
obtained from the noise map for snapshot i. This com- 
bination optimizes the signal to noise ratio of the final 
image. Beam patterns are calculated independently for 
each 7.68 MHz channel. The same weighting scheme was 
used to combine the "clean components" and residual 
maps of the individual snapshots. 

A v ariant of the Hogbom CLEAN algorithm (jHogbomI 
Il974f ) was used to further deconvolve the integrated 
residual maps, using a position-dependent point spread 
function calculated using the same weighting scheme: 
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(7) 

where Pint.j is the point spread function for a source at 
position {9j ,(j)j), Bi is the primary beam pattern for the 
ith snapshot with a PSF given by Pi, and the peak of 
the function is normalized to unity. CLEAN components 
were selected by choosing the pixel in the residual map 
with the largest signal-to-noise ratio (SNR, determined 
by dividing the residual map by its noise map). The PSF 
was scaled to a peak of 10% of the flux of the pixel. The 
images were restored with a Gaussian beam determined 
by a fit to the weighted average of the individual snapshot 
point spread functions at the field center. 

6.5. Averaging 

In order to increase the signal to noise ratio and im- 
age fidelity for source detection and characterization, 
the individual 7.68 MHz maps, after deconvolution and 
restoration, were averaged together. An approximate 
flux scale for the maps was first set by scaling the surface 
brightness of the maximum pixel at the location of Hy- 
dra A to a value of 296 x {iy/150 MHz)-°-9i Jy beam"^ 
(this model was derived from fitting a power law to other 
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low-frequency measurements). For each field, 30.72 MHz 
bandwidth maps centered at 123.52 MHz, 154.24 MHz, 
and 184.96 MHz were each made from four 7.68-MHz 
maps. Before averaging, the 7.68-MHz map fluxes were 
scaled to the averaged map frequency using an assumed 
spectral index of a = —0.8 (where S (xv"). The averages 
were computed in a weighted sense using the integrated 
primary beam weights from each map. A full-band (92.16 
MHz bandwidth) weighted average map was made from 
the three 30.72-MHz bandwidth maps after scaling them 
to a common reference frequency of 154.24 MHz, again 
using a spectral index of a = —0.8. The portion of each 
field within 25° of the field center was used for the sub- 
sequent analysis. 

6.6. Source Extraction 

Sources were identified in each full-band, i.e., 92.16 
MHz wide, averaged map using an automated source ex- 
traction pipeline. The first step of this pipeline was the 
calculation of a position-dependent noise map using the 
method described in Section [6^ The full-band map was 
then divided by the noise map to produce a SNR map. 
An iterative process then was initiated by identifying 
contiguous regions of pixels above a certain SNR detec- 
tion threshold and, for each such region, defining a fitting 
region that extended beyond the set of connected pixels 
by several synthesized beamwidths. A two-dimensional 
Gaussian fit was performed on the corresponding region 
in the full-band map. The parameters determined in 
each fit included the background level, peak position, 
peak amplitude, major axis width, minor axis width, 
and position angle. If the fit converged to a Gaussian 
centered within the fitting region, then the source was 
subtracted from the map, and the extracted source pa- 
rameters were recorded. After fitting all regions identi- 
fied for a certain SNR level, the detection threshold was 
reduced and the process was repeated. Regions above the 
detection threshold for which a fit failed to converge are 
refit in subsequent iterations at lower detection thresh- 
olds (where the fitting regions are typically larger in size). 
For completeness, sources were extracted down to a de- 
tection SNR threshold of 3. It should be noted that in 
this fitting procedure, each region that is fit by a sin- 
gle two-dimensional Gaussian is taken to correspond to 
a separate source. Sources that are too close together to 
be resolved into separate components will be fit with a 
single component and erroneously taken to be a single 
source that is a "blend" of the two components. 

Sources that were identified in the full-band average 
map were then extracted from each of the 30.72-MHz 
bandwidth maps. The sources were sorted by their detec- 
tion signal-to-noise level and, for each of the three maps, 
were fitted in descending order. For each source, the po- 
sition and shape (axial ratio and position angle) param- 
eters of the Gaussian fitting function were held fixed to 
the values determined in the full-band map extraction, 
while the peak value, background level, and a scaling 
factor for the widths of both the major and minor axes 
of the Gaussian were allowed to vary. This procedure 
was performed for all sources. When the fit successfully 
converged, the best-fit model was subtracted from the 
sub-band map. A total of 908 sources were extracted in 
the Hydra A field and 1100 sources were extracted from 
the EoR2 field. 



6.7. Astrometric Corrections & Flux Calibration 

We compared the positions and relative fiuxes of the 
sources identified in the full-band maps with the posi- 
tions and fluxes of possible counterparts in other cata- 
logs. These comparisons formed the basis for astromet- 
ric corrections, and for the determination of the over- 
aU MWA flux scale. Counterparts to MWA-32T sources 
were identified at 408 MHz by locating sources in the 
MRC which were within 15' of an MWA source. Al- 
though we expect our astrometric accuracy to be much 
better than 15', this value was chosen to be compara- 
ble to the size of the MWA-32T synthesized beam major 
axis full width at half maximum response, viz., 13' for 
the Hydra A field and 14' for the EoR2 field, and to 
allow for some degree of systematic error in the MWA 
source positions. To avoid possible blending issues, we 
only considered an identification to be secure when there 
was precisely one source in the MRC within 30' of the 
(pre-adjustment) MWA source position. A total of 419 
sources were matched uniquely to the MRC in the Hy- 
dra A field and 520 sources were matched in the EoR2 
field. 

An astrometric correction was then calculated by al- 
lowing for a linear transformation of the MWA source 
coordinates in order to minimize the positional differ- 
ences between corresponding MWA and MRC sources. 
The transformation permits offsets in both right ascen- 
sion and declination as well as rotation and shear with 
respect to the field center. Optimal transformation pa- 
rameters were determined by performing weighted least- 
squares fits. The results from initial fits indicated that 
there were errors in the (Earth-referenced) coordinates 
of the MWA tiles and in the conversion of coordinates 
in the maps from the epoch of observation frame to the 
J2000 frame. These errors were then corrected. The fi- 
nal fits were found to be consistent solely with offsets of 
the MWA source coordinates, with no shear or rotational 
effects. We therefore applied offsets of Aa = —0.6' and 
AS = 1.6' to the positions of the sources in the Hydra A 
field, and of Aa = 2.2' and Ad = 1.9' to those in the 
EoR2 field. We believe the coordinate offsets are likely 
due to a combination of the effects of structure in Hy- 
dra A and to ionospheric refraction. Figure [2] shows the 
post-offset-correction differences in position between cor- 
responding MWA and MRC sources in the EoR2 field. 
Histograms of residual positional differences from both 
the Hydra A and EoR2 fields are plotted in Figure [31 For 
this figure the difference for each source is normalized by 
the expected error based on the signal to noise ratio for 
the intensity and the 32T synthesized b eamwidth (with 
the assumption of a circular source, see ICondonI [19971. 
The residual position differences are generally consistent 
with those expected, even though there are a small num- 
ber of large position differences. 

A final flux scale was set for each map using the MWA 
sources which had counterparts in both the MRC and 
Culgoora source lists. Only sources in the MWA catalog 
above a detection SNR threshold of 5 were used. Hy- 
dra A was excluded for reasons discussed in Section 17.11 
A prediction for each source was obtained by fitting a 
power law spectrum to the 408 MHz MRC flux and the 
80 MHz and 160 MHz Culgoora fluxes. Under these cri- 
teria, measurements in all 3 bands were found for a total 
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Fig. 2. — Spatial oflFsets between the positions of MWA sources 
in the EoR2 field and matched sources in the Molonglo Reference 
Catalog (Large et al. 1981^ An overall coordinate system shift has 
been removed. 
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Fig. 3. — Histograms of the normalized right ascension, a, (a) and 
declination, 5, (b) errors for the extracted MWA sourc es relative to 
matc hed sources in the Molonglo Reference Catalog l ILarge et al.l 
119811 1. The standard deviations, ao:, and ag, are calculated fol- 
lowing [Condo^ l|1997l ) with the simplifying assumption of circular 
source geometry Assuming Gaussian error properties, the residual 
distribution should approximate a standard normal distribution, 
which is over-plotted with a dashed line. Hydra A is omitted from 
the histograms. 

of 64 uniquely matching sources in the Hydra A field and 
81 uniquely matching sources in the EoR2 field. Using 
these flux predictions, a flux scale correction was calcu- 
lated of the form: 



•^calibrated — ^u ^ *^uncalibratcd- 



(8) 



These calibration terms were calculated independently 
for each of the three sub-band maps as well as for each av- 
eraged map for each field. For the Hydra A field, the cal- 
culation yeilded Ci, ~ 1.26 for the full-band average map, 
1.17 for the 123.52 MHz map, 1.20 for the 154.24 MHz 
map and 1.18 for the 184.96 MHz map, and for the EoR2 
field, the calculation yielded Ci, ~ 1.24 for the full-band 
average map, 1.19 for the 123.52 MHz map, 1.23 for the 
154.24 MHz map and 1.17 for the 184.96 MHz map. The 
residuals after applying these flux scale corrections were 
analyzed to determine if additional biases were present as 
a function of position in the image (biases would poten- 
tially be seen, e.g., if the assumed primary beam model 
was incorrect). An example of these plots is shown in 
Figure |4l which shows no evidence for a radially increas- 
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Fig. 4. — The fractional differences between predicted and mea- 
surc d flu xes for the EoR2 field at 154.24 MHz (described in Sec- 
tion [BTTll, where AS is defined as the MWA-32T measured fiux 
minus the flux predicted from fitting MRC and Culgoora measure- 
ments. The error bars are derived from combining the RMS noise 
in the MWA map at the source position in quadrature with the flux 
prediction uncertainty. The differences are displayed as a function 
of distance from the field center, in order to assess the presence of 
any radial biases in the MWA flux measurements. No significant 
bias was found out to large distances from the field center. 

ing flux bias. 

The magnitude of the post-correction residual differ- 
ences between the MWA measurements and the pre- 
dicted fluxes for calibration sources are still larger on 
average than expected under the assumption that the un- 
certainty in each MWA-32T flux measurement is due to 
the RMS map noise at the source location, and that the 
uncertainty in the predicted flux of each source is prop- 
agated for the power-law fitting procedure. This excess 
in the differences could be due to errors in the spectral 
model for the calibration sources, temporal variability of 
the sources, or to as yet unidentified errors. Since this 
is the first work based on MWA-32T data to report the 
fluxes of a large number of sources, we make the conser- 
vative assumption that these excessively large residuals 
are due solely to errors in the MWA-32T flux measure- 
ments. We assume the flux uncertainties follow a Gaus- 
sian distribution which includes the effects of the RMS 
map noise added in quadrature with an additional com- 
ponent proportional to the measured flux of the source: 
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where ctmwa is the la flux uncertainty for a particular 
source, p is the fractional flux uncertainty, S'mwa is the 
measured source flux, and CMap is the RMS map noise 
at the position of the source. We evaluated the standard 
deviation, ajj, of the fractional flux difference, D, where 
D is calculated as: 
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We then solved for the fractional uncertainty in the MWA 
measurements which would be needed to account for the 
magnitude of the measured value of ao '■ 
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We calculated the average value of /? separately for each 
held. For the higher- frequency maps, we found that the 
values of j3 were much larger far from the field center 
where the primary beam approaches the first null; for 
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these maps the sources were separated into inner and 
outer region sets using a cutoff of 18°, and /3 was calcu- 
lated separately for each region. Using these results, we 
assign fractional flux uncertainties of 30% for all sources 
in the full-band average maps, 35% for all sources in the 
123.52 MHz maps, 35% for sources in the inner region 
of the 154.24 MHz maps, 60% for sources in the outer 
region of the 154.24 MHz maps, 35% for sources in the 
inner region of the 184.96 MHz maps and 80% for sources 
in the outer region of the 184.96 MHz maps. These frac- 
tional uncertainty values are applied to all sources in the 
catalog by adding them in quadrature to the map RMS 
values as described in Equation [Sj 



expression takes the form of a broken power la'wrn: 



7. RESULTS 
7.1. Radio Maps and Source Catalog 

The full band average maps of the Hydra A and EoR2 
fields are displayed in Figure [S] and Figure HI These 
images overlap partially. Together, they cover ~2700 
square degrees. The synthesized beam for the Hydra A 
field has a major-axis width of 19' in the 124.52 MHz 
map, 14' in the 154.24 MHz map, 12' in the 184.96 MHz 
map, and 13' in the full band average map. For the 
EoR2 field, the major-axis beam widths are 18' in the 
124.52 MHz map, 16' in the 154.24 MHz map, 13' in the 
184.96 MHz map, and 14' in the full band average map. 

The bright radio galaxy Hydra A is the dominant 
source in these maps. The flux of the source is mea- 
sured to be 710 ± 210 Jy in the full band average map 
of the Hydra A field, and 550 ± 170 Jy in the full band 
average map of the EoR2 field. These measurements of 
Hydra A are significantly brighter than expected based 
on previous measurements: the Culgoora 160 MHz mea- 
surements give a fiux of 243 Jy, and a prediction based off 
of the Culgoora and MRC measurements (as described in 
Section 16. 7|) gives a fiux of 284 Jy. Although Hydra A is 
slightly extended in the MWA-32T maps, the structures 
seen in th e previous low freq uency maps of Hydra A pre- 
sented by iLane et all ()2004[ ) are below the scale of the 
MWA synthesized beam. We note that the MWA-32T 
array has a significantly more compact uv distribution 
than the VLA, Molonglo or Culgoora telescopes. A flux 
from Hydra A above what is expected from the Culgo- 
ora measurements is also note d in measurements with 
the PAPER array pacobsll2011l private communication) , 
which has a similar baseline distribution to that of the 
MWA-32T. Consequently we did not include Hydra A 
in the final fiux scale calibration procedure described in 
Section 16.71 

The behavior of the fitted RMS noise in these images 
is illustrated in Figure [3 which shows annular averages 
as a function of distance from the field center. We can 
estimate a lower limit to the RMS noise in each map by 
calculating the classical source confusion limit. Note that 
this differs fr om the sidelobe confusion limit; see, e.g. 
lCondonlll974l for a rigorous disc ussion of classical source 
confusion in radio telescopes. iDi Matteo et al.l (|2002[ ) 
present a model of the r adio source counts derived from 
the 6C (|Hales et al.l [l98& ) catalog at 150 MHz. Their 



dn (4000 {T%y^'^'^ sources Jy"^ sr"\ 



dS 



4000 (t%)-°''' (tt?)"'"' sources Jy-i 



S > So 
S <Sa 
(12) 

where Sq — 0.88 Jy. Integration of this expression gives 
the source density in each beam above a minimum flux 
value, ^min: 



PS 



41n2 



dn 
dS 



dS sources beam 



(13) 



where ps is the source density in units of sources per 
beam, and 9 is the FWHM synthesized beam size. The 
source confusion limit then corresponds to the flux, S'min, 
for which the source density approaches one source per 
synthesized beam area (typically maps are considered to 
be source confusion limited when they have a source den- 
sity of greater than 1 source per 10 synthesized beams). 
The average RMS noise (Figure [7]) reaches a minimum 
value of ^160 mJy for the full band average map of the 
EoR2 field. Using the corresponding fiux threshold for 
an unresolved source of S'min — 160 mJy in combination 
with the EoR2 full band map synthesized beam area gives 
PS ^ 0.30 sources per synthesized beam, or, in other 
words, one source of at least 160 mJy in roughly every 
three synthesized beams. Estimates for the other average 
and sub-band maps give expected source densities of one 
source per every 3 to 6 synthesized beams . Thus, source 
confusion is likely the limiting source of noise in the cen- 
tral region of these maps. This explains the relatively flat 
nature of the central noise floor seen in Figure [71 Near 
the edges of the images, not far from the flrst null of the 
primary beam, we expect the noise to be dominated by 
receiver noise, and to scale with the inverse of the pri- 
mary beam power pattern. This is also seen in Figure [7] 
(the frequency dependence of the noise curves illustrates 
the frequency dependence of the primary beam). 

A catalog was constructed from the 2008 detections at 
a detection SNR threshold > 3 of potential sources in the 
EoR2 and Hydra A fields. In the cases where there were 
detections at corresponding celestial positions in the two 
fields, the measurement where the source was closer to 
the observation field center was retained, resulting in a 
list of 1526 unique source detections. The quality of this 
source list is assessed in Section 17.21 and Section 17.31 as 
a function of the detection SNR level. The 655 sources 
detected at an SNR level >5 in the detection images are 
reported in TableU The 871 sources with 3<SNR<5 are 
considered to be less reliable detections. A list of these 
candidate detections can be obtained by contacting the 
authors. 

7.2. Reliability of the Source List 

The reliability of the identified sources was evaluated 
through comparison with the fiux limited sample from 



^^ The notation of lDi Matteo et al.l IJ2002I ) is ambiguou s, and is 
interp reted differently by several authors. We note that Li dz et al.l 
II2008I ) quote a modified form of the expression, which affects the 
normalization of the power law. Although the Lidz ct al. (2008) ex- 
pression fits the 6C source counts slightly better at high flux values, 
the formulation presented in this paper fits the data better through- 
out th e entire flux range. The expression from Di Matteo ct al] 
1)20021 ) has been corrected with an additional minus sign to make 
it continuous across the transition at S = So- 
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Fig. 5. — Full bandwidth synthesis map of the Hydra A field, with a 50° 
in Section |6] and served as a basis for source identification in this field. 



diameter. This map was produced using the pipeline described 



the MRC, the VLSS, and maps from the TGSS. The 
MRC source Ust has a well-defined completeness flux 
limit and covers our entire field; however it gives fluxes 
at a different frequency (408 MHz) and does not go quite 
as deep as our survey. The VLSS covers a portion of our 
fields at a lower frequency than the MWA (74 MHz), and 
provides a useful complementary assessment. The TGSS 
maps are at the same frequency as the MWA observa- 
tions (150 MHz), but the maps which have been released 
to date only cover a small fraction of our fields and are 
based on a significantly different sampling of the visibil- 
ity function due to the different array baselines. 

In order to assess the MWA-32T catalog reliability, we 
first evaluate the detectability of an MWA source in the 
external comparison survey. The MWA full band aver- 
age flux is extrapolated to the relevant frequency using a 
spectral index of a = —0.8, and the extrapolated value is 
then compared to the parameters of the comparison map 



or catalog to evaluate whether it meets the detection cri- 
teria for that survey. If the source is deemed detectable, 
then we search for a companion source in that catalog 
or map to see if the source was actually detected by the 
other survey. Under the assumption that the other sur- 
veys are complete, this allows us to assess how many 
spurious sources are present in the MWA catalog. We 
define the reliability as: 



R 



-^detected/A^detectable — / 



(14) 



where R is the fraction of MWA sources we believe to be 
reliable, A'dctcctabic is the number of MWA sources which 
we believe should have been detectable in the compari- 
son survey, A'^dotocted is the number of detectable MWA 
sources for which we found counterparts in the other sur- 
vey, and / is the false source coincidence fraction. We 
determine / by calculating the source density of the com- 
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Fig. 6. — Full bandwidth synthesis maps of the EoR2 field, with a 50° 
Section |6] and served as a basis for source identification in this field. 

parison survey in the MWA fields, and use our counter- 
part matching criteria to estimate the probability that 
a randomly chosen sky location will lead to an associa- 
tion with a source in the comparison survey. This anal- 
ysis was performed for different MWA source detection 
thresholds. 

For the reliability compariso n with the MRC, w e used 
the completeness limit of 1 Jy (jLarge et al.lll98l[) to as- 
sess the detectability of the extrapolated MWA source 
fluxes. An MRC counterpart is associated with the MWA 
source if it is within 10' of the MWA source position. 
Based off of the counterpart search radius and source 
densities in the MRC field, we estimate a false coinci- 
dence chance of 4%. A fixed fiux completeness limit is 
not given for the VLSS, however Cohen ct al. (2007) note 
that for a typical VLSS RMS of 0.1 Jy beam-\ the 50% 
point source detection limit is approximately 0.7 Jy. We 
then assume the VLSS is complete down to a flux level 



diameter. This map was produced using the pipeline described in 

of 1 Jy, and we again use a 10' source association radius 
in the reliability calculation. At the present time, the 
VLSS catalog does not cover the entire combined EoR2 
and Hydra A region that we have surveyed. To ensure 
we are only including sources in the VLSS survey area, 
we only analyzed MWA sources above a declination of 
5 = —25°. We estimate a false coincidence chance of 
18% for the VLSS matching. 

The cumulative and differential catalog reliabilities are 
listed as functions of the MWA detection level in Ta- 
ble [3] and Table U) We view these reliability estimates 
as lower limits, particularly at the lower flux levels, be- 
cause our assessments of comparison survey detectability 
do not take into account errors in the source flux extrap- 
olation or source time variability. MWA sources which 
are erroneously calculated as detectable will not, in gen- 
eral, lead to detections of counterparts in the comparison 
catalog, whereas sources erroneously calculated as unde- 
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TABLE 2 
Detected sources in the the Hydra A and EOR 2 fields 



Name 


RA 


DEC 


Savg 


•S'l23.52 


5l54.24 


S'l84.96 


Field 


rpc 


Detection 






















SNR Level 


J0747- 


-1854 


07h47m05s 


-18°54'12" 


8.8±2.8 


10.1±3.7 


9.2±6.0 


8.9±8.5 


HydA 


22.9° 


8.6 


J0747- 


-1919 


07h47m27s 


-19°19'33" 


22.2±6.7 


25.2±8.9 


23.8±14.5 


15.1±13.0 


HydA 


23.0° 


22.4 


J0751- 


-1919 


07h51m20s 


-19°19'30" 


7.1±2.2 


8.1±3.0 


8.5±5.3 


4.7±5.1 


HydA 


22.1° 


8.7 


J0752- 


-2204 


07h52m30s 


-22°04'38" 


4.4±1.5 


4.7±1.8 


7.1±4.5 


14.4±12.1 


HydA 


22.7° 


5.6 


J0752- 


-2627 


07h52m30s 


-26°27'43" 


8.4±2.7 


9.8±3.7 


3.6±3.1 




HydA 


24.7° 


7.0 


J0757- 


-1137 


07h57ml0s 


-11°37'10" 


3.9±1.3 


4.7±1.8 


5.1±3.3 


2.5±2.5 


HydA 


19.8° 


6.0 


J0802- 


-0915 


08h02ml8s 


-09°15'32" 


3.4±1.2 


5.8±2.2 


2.1±1.5 


1.9±1.8 


HydA 


18.8° 


5.0 


J0802- 


-0958 


08h02m34s 


-09°58'55" 


8.9±2.8 


12.4±4.4 


6.0±3.7 


9.1±7.3 


HydA 


18.6° 


12.6 


J0803- 


-0804 


08h03m60s 


-08°04'48" 


4.5±1.4 


6.5±2.4 


2.4±1.6 


6.7±5.4 


HydA 


18.7° 


7.0 


J0804- 


-1244 


08h04ml7s 


-12°44'32" 


4.6±1.5 


4.8±1.8 


5.2±3.2 


4.5±3.7 


HydA 


18.0° 


9.3 


J0804- 


-1726 


08h04m42s 


-17°26'41" 


4.9±1.5 


5.8±2.1 


4.5±2.8 


3.1±2.6 


HydA 


18.5° 


9.1 


J0804- 


-1502 


08h04m53s 


-15°02'54" 


2.8±0.9 


2.2±1.0 


4.6±2.9 


l.Oitl.l 


HydA 


18.0° 


6.4 


J0805- 


-0100 


08h05m30s 


-01°00'12" 


9.4±2.9 


10.8±3.9 


8.3±5.2 


6.4±5.4 


HydA 


21.1° 


10.4 


J0805- 


-0739 


08h05m40s 


-07°39'22" 


3.5±1.2 


4.4±1.7 


3.7±2.4 


1.9±1.9 


HydA 


18.4° 


5.5 


J0806- 


-2204 


08h06m26s 


-22°04'43" 


3.7±1.2 


4.2±1.6 


3.8±2.4 


3.6±3.1 


HydA 


19.8° 


6.2 



Note. — The flux of each source detected in the MWA full-band averaged maps is presented along with the flux 
measured in each 30.72 MHz sub-band. Duplicate sources in the region of overlap of the two flelds are not listed. Missing 
data indicates that the automatic source measurement algorithm failed to converge in a flux fit for that source in the 
sub-band map. The field from which each source measurement comes from is listed, along with the distance of the source 
from the center of the field (rpc)- We expect systematic errors to be larger for sources far from the field center. The 
"Detection SNR Level" indicates the signal-to-noise ratio at which the source was detected in the full-band averaged map. 
Section 17. 21 discusses the reliability of the catalog at different detection SNR levels. This list includes sources identified above a 
detection SNR threshold of 5. The full source list of all sources above a detection SNR threshold of 3 is available from the authors. 

Table \2\ is published in its entirety in the electronic edition of the Astrophysical Journal. A sample is shown here to 
illustrate its content. 
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Fig. 7. — Radial dependence of the calculated RMS noise in the 
Hydra A (a) and EoR2 (b) field images. The RMS is plotted for the 
three sub- band maps, with central frequencies of 123.52 MHz (solid 
Unes), 154.24 MHz (dashed lines), and 184.96 MHz (dotted line) as 
well as for the full- band averaged maps (dot-dashed lines). The val- 
ues shown are the medians of 1° wide radial annuli. The points are 
connected with lines for clarity. The frequency dependent primary 
beam shape is evident from the increasing RMS at large distance 
from the field center, and the confusion limited nature of the maps 
is illustrated by the fiat RMS profile near the center of the field. 
The minimum RMS approaches a value of ~130 mjy/beam. 

tectable will be omitted from the analysis and, therefore, 
will not be included in the calculation of the reliability 
ratio. These catalog comparisons imply a reliability of 
> 99% for sources detected above a detection SNR of 5. 
For the reliability comparison using the TGSS, we used 
the maps from TGSS Data Release 2 available at the time 
of our analysis. Although the baseline distribution of the 
GMRT is substantially different from that of the MWA- 
32T, the GMRT has a compact central array consisting 
of fourteen antennas within an area of radius 500 meters 
(jSwarup et al.l 119911) that leads to substantial overlap 
with the MWA-32T regarding the region of the uv plane 
that was sampled. Twenty-seven TGSS fields overlapped 
the EoR2 field; none overlapped the Hydra A field. We 



modeled the effects of source blending by the large MWA 
beam by convolving the CLEANed and restored TGSS 
maps to a Gaussian full width at half maximum res- 
olution of 12'. Care was taken in the convolution to 
preserve the flux density scale. The resulting convolved 
images have typical RMS surface brightness fluctuations 
of '^0.4 Jy/beam in regions free of sources. Sources in 
these maps were identified by taking all pixels above 4cr 
and associating a source with each island of bright pix- 
els. The resulting TGSS source catalog was compared to 
those sources in the MWA list with flux density greater 
than the 4cr level in the TGSS field, and thus expected 
to be detected in the TGSS field. Pairs of sources in 
the two catalogs coincident within 10' were recorded as 
sources detected in both surveys. MWA sources with- 
out a TGSS counterpart were counted as non-detections. 
Reliability values are presented in Tables |3] and S) These 
values are ratios for each MWA detection SNR bin of 
the number of TGSS detections to the number of MWA 
sources expected to be detected in the TGSS field, and 
are, of course, a function of the threshold chosen in the 
TGSS maps. As an example of the TGSS comparison we 
present Figure [HI which plots the positions of MWA and 
MRC sources on a grey scale image of a convolved TGSS 
field. It is important to note that these reliability esti- 
mates solely test for the presence of a source coincident 
with the reported position, and do not speak to the fi- 
delity of the fluxes of these sources or whether the MWA 
sources are due to single objects or blends of multiple 
fainter objects. 

7.3. Completeness of the Source Catalog 

We used the Culgoora source list to assess the com- 
pleteness of the MWA catalog presented in this paper. 
We chose the Culgoora list because it includes observa- 
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TABLE 3 

Cumulative Source Reliability 



Detection SNR 
Threshold 



MRC Reliability 

-'* detected /-'* detectable 



R 



VLSS Reliability 

-^ 'detected/-' * detectable 



R 



TGSS Reliability 

""^ detected/-^ 'detectable 



> 3 


488/589 


82% 


1215/1312 


91% 


183/197 


93% 


> 4 


421/444 


95% 


826/839 


98% 


132/133 


99% 


>5 


349/357 


98% 


575/579 


99% 


85/85 


100% 


> 7 


257/259 


99% 


325/326 


100% 


49/49 


100% 


> 10 


167/167 


100% 


173/173 


100% 


26/26 


100% 



Note. — The reliability is assessed by comparing the MWA source list with the MRC catalog, VLSS catalog, and 
convolved TGSS maps as descriged in Section |7.2| The MRC and VLSS comparisons are made by extrapolating the MWA 
source flux to 408 MHz and 74 MHz respectively, assuming a spectral index of a — —0.8, to assess the dctcctability of 
the MWA source in the catalogs. The TGSS results are based on using sources in convolved TGSS maps above 4(7 
significance. The reliability percentages, R, have been corrected for false positives. 



TABLE 4 

Differential Source Reliability 



Detection SNR 


MRC Diff. Reliability 


VLSS Diff. 


Reliability 


TGSS Diff. Reliability 


Threshold (DT) 


-'Vdctcctcd/-'VdctcctabIc R 


■'Vdotcctcd/-'Vdct< 


ctabic 


R 


-'Vdctoctod /'^detectable 


R 


3 < DT < 4 


67/145 44% 


389/473 




78% 


51/64 


80% 


4 < DT < 5 


72/87 82% 


251/260 




96% 


47/48 


98% 


5 < DT < 6 


43/48 89% 


142/144 




98% 


20/20 


100% 


6 < DT < 7 


49/50 98% 


108/109 




99% 


16/16 


100% 


7 < DT < 8 


33/34 97% 


68/69 




98% 


15/15 


100% 


8 < DT < 9 


31/32 97% 


44/44 




100% 


7/7 


100% 


9 < DT < 10 


26/26 100% 


40/40 




100% 


1/1 


100% 



Note. — The reliability is assessed by comparing tlie MWA sourec list with the MRC catalog, VLSS catalog, and 
convolved TGSS maps for various ranges of detection SNR threshold (DT) as discussed in Section 17.21 The MRC and 
VLSS comparisons are made by extrapolating the MWA source flux to 408 MHz and 74 MHz respectively, assuming a 
spectral index of a — —0.8, to assess the dctcctability of the MWA source in the catalogs. The TGSS results are based 
on using sources in convolved TGSS maps above 4(T significance. The reliability percentages, R, have been corrected for 
false positives. 



tions done at 160 MHz, a frequency not far from the 
midpoint of the MWA band and because its synthesized 
beam is similar in size to that of the MWA-32T. For each 
Culgoora source within a field observed by the MWA 
we used the position-dependent RMS noise in the MWA 
full-band averaged maps to evaluate its detectability in 
the MWA map. Culgoora sources which should be de- 
tectable above a specified SNR level in the MWA images 
were then checked for a matching source within 10' in 
the MWA catalog. Since we used the Culgoora source 
list to assess the completeness, the results are only valid 
down to a level comparable to the lowest Culgoora fluxes 
of ~ 1.2 Jy. The completeness ratio was calculated sim- 
ilarly to the reliability described in Section [721 



C = 



m 



-i/N, 



detected/-^ * detectable 



-/ 



1-/ 



(15) 



where C is the completeness percentage, A^dctectabie is 
the number of Culgoora sources which we believe should 
have been detectable in the MWA source list , A'^detected 
is the number of detectable Culgoora sources for which 
we found counterparts in the MWA list, and / is the false 
source coincidence fraction calculated from the MWA 
catalog source density using the 10' source matching cri- 
teria. The results are presented in Table [S] 

The completeness was analyzed separately for sources 
within inner and outer regions separated by a circle of 



radius r — 18° around the field center. All Culgoora 
sources within the inner region which did not have a cor- 
responding detection in the MWA source list were in- 
spected and found to coincide with a local maximum in 
the map, implying the completeness is limited by the ro- 
bustness of the source extraction algorithm and the flux 
calibration rather than the intrinsic map quality. It is 
important to note that because the MWA maps have a 
sensitivity that varies strongly across the field, this com- 
pleteness value does not specify a flux limit to the cata- 
log, but rather assesses the efficacy of the source extrac- 
tion. As with the above reliability estimate, variability 
and incompleteness act to make this a lower limit on the 
true completeness. 

7.4. Source Counts and Correlation Function 

Radio source counts provide another useful diagnos- 
tic test to assess the quality of the catalog and co nsis- 
tency with previ o us wo rks. As discussed in Section 17.11 
iDi Mattco et all (I2002D fi t the 151 MHz 6C survey re- 
sults of Hales et alj ( 1988D to obtain a power law model 
for radio source counts. The fit is valid up to ^ 10 Jy, 
but the actual 6C counts fall somewhat below the fit at 
the high end of the range. Using the MWA-32T catalog 
generated from the EoR2 field, we calculate the differ- 
ential source counts, using the noise map of the field to 
correct for the effects of sensitivity variations on the ef- 
fective survey area for different flux values (i.e., bright 
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TABLE 5 
Source List Completeness 



Field r < 18 
Name A^dotcctcd/A^dctcctablc 




Completeness 


■'Vdctcc 


■r > 18° 

tcd/-'Vdctcctablc 


Completeness 


Detection SNR Level > 5 

Hydra A 56/63 

EoR2 72/77 

Detection SNR Level > 3 

Hydra A 62/63 

EoR2 75/77 






89% 
93% 

98% 
97% 




36/44 
49/58 

52/67 
61/66 






82% 
84% 

77% 
92% 


Note. — The completeness 
I1995D. as described in Section 


as assessed by a comparison wi 
17.31 The minimum source flux 


th the Cu 
in the Cr 


goora 
Igoora 


source list dSleel 
list is ~1.2 Jy, 



SO these results arc only valid for sources brighter than this level. We view these completeness 
estimates as a lower limit on the catalog completeness — source variability or flux errors in 
the Culgoora measurement will decrease the calculated completeness ratio. Due to the varying 
sensitivity across the MWA fleld, the completeness is calculated relative to the local noise in the 
MWA source detection map, rather than an absolute flux level. The completeness percentages 
have been corrected for false positives as described in Section |7.3l Analyzing the source counts 
in the fleld (sec Section |7.4| l indicates that the source list is complete above ~2 Jy. 
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Fig. 8.— Grey scale image of a TGSS field (field R33D18), 
with positions of MWA and MRC sources overlaid. Only pix- 
els with SNR > 4 are plotted, and the mapping of pixels to 
grey scale is shown by the scale in Jy beam" at the top of the 
image. Positions of MWA sources with Detection Threshold > 
5 are plotted with a square, of MWA sources with 4 < 
Detection Threshold < 5 are plotted with a circle, of MWA source 
with 3 < Detection Threshold < 4 with a triangle, and of MRC 
sources with an X. All MWA and MRC sources in this field have a 
counterpart in the TGSS image. There are five sources in the field 
that are detected in the MWA and TGSS surveys, but not in the 
MRC. 



sources can be detected over a larger area than faint 
sources). No Eddington bias correction is applied (to 
correct for the artificial enhancement of faint sources 
due to noise in the map) and the error bars are calcu- 
lated from the square root of the number of counts in 
each bin. These source counts are shown in Figure [SI 
alo ng with the expected so urce counts from integrating 
the i Di Matteo et al.l (|2002[ ) model. We note that more 
sophisticated models of 150 MHz source counts mod- 



els have been desc ribed by, e.g., IWilman et al.l ()2008f ) 
and lJacksonI (J2005D . These models have different behav- 
ior at the sub-Jansky flux levels that have been probed 
by high resolution, deep, nar rower field-of -view stud- 
ies such as those descr i bed in llntema et al] {2011) and 
llshwara-Chandra et al.l (I2010D. We compare our results 
with the IDi Matteo et al.l ()2002l ) source counts model be- 
cause it is commonly used as a basis for studies of EoR 
foregrounds and sensitivities. 

The results fr om the catalog p r esente d in this work 
agree with the IDi Matteo et al.l ()2002D model above 
flux levels of ~2 Jy Below this level, the MWA- 
32T source counts diverge from the model, likely be- 
cause of the incompleteness of the MWA source ex- 
traction for low flux sources. A power law fit to the 
EoR2 field source counts above 2 Jy yields dn/dS = 
(3500 ± 500)(S'/1 jy)-2.59±o.09 sources Jy"^ sr"! for 
sources with a detection SNR greater than 5. Fit- 
ting for a power law to the full list of sources in the 
field down to a detection SNR of 3 yields dn/dS = 
(5700 ±700) (5/1 jy)-2.76±o.08 sources Jy"^ sr"!. 

As an additional test for systematic effects, we have 
constructed the angular two-point correlation function, 
w{9), of the sources in our catalog. This correlation func- 
tion can show systematic effects that manifest themselves 
on char acteristic angular scales in the catalog ~ see, e.g. 
IBlake fc Wall (20021. We measur e w{6) using the esti- 
mator defined bv iHamiltonl (|1993[ ): 



.,,.iMM_, 



(16) 



where DD{9) is the measured angular autocorrelation 
function from the MWA source catalog, RR{9) is the 
autocorrelation function calculated using a simulated 
"mock" catalog, and DR{9) is the cross-correlation be- 
tween the MWA and the mock catalog. We gener- 
ate an ensemble of 100 mock catalogs and evaluate 
the correlation function with each one separately in 
order to produce a set of normally distributed esti- 
mates of w{9). Each mock catalog is produced us- 
ing an approach developed to simulate point sources 
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Fig. 9. — DiflFerential source counts histograms from the MWA EoR2 field, calculated using both the high reliability catalog (a) and 
the full list of source candidates (b). The noise maps were used to correct for the effective area surveyed in each bin. Poisson error 
bars are assigned based on the number of counts in each bin. No Eddington bias correction is app lied. The source counts model from 
IDi Matteo et al.i ^2002) is shown for reference with a dashed line. We believe the deviations from the lDi Matteo et al.l IJ2002I ) model below 
5 ~ 3 Jy are not due to an intrinsic change in the source counts distribution at this scale, but are instead due to incompleteness of the 
source catalog in this flux range. 



at CMB and FIR frequencies ('Argiicso et al.1 120031: 
I Gonzalez-Nuevo ~et al . 2005), but tailored spec ifically fo r 
the MWA experiment (de O liveira-Costa et aJllin prep.f ). 
i.e., we drew sources from the observed MWA-32T 
source counts distribution described above in accor- 
dance with the expected low-frequency source clus- 
tering statistics Ide Oliveir a-Costa fc Capodilupol 120101 : 
Ide Oliveira-Costa fc Lazioll2010D . On the angular scales 
probed by this survey, no observable clustering is ex- 
pected. By constructing the mock catalogs in this man- 
ner, and correlating them with the observed distribution, 
the resulting estimate of w{9) identifies any unexpected 
correlation which may be due to systematic errors in our 
survey or in our catalog construction procedure. Fig- 
ure 1 101 shows our measurement o{w{9) above a flux limit 
of S" « 3 Jy (black squares). Distances between the ob- 
served and/or simulated sources are measured in bins 
of ^ 1°, which is substantially above the MWA reso- 
lution. The mean value of w{9) is shown, along with 
uncertainties derived from calculating the covariance be- 
tween w{6) bins in the mock catalogs. As expected, 'w{9) 
is consistent with zero, implying that there is no excess 
correlation in our catalog. 
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Fig. 10. — Measured angular correlation function, w{9), of the 
full MWA catalog. The angular correlations are calculated using 
sources with fluxes above 5^3 Jy. The top panel shows w{8) 
calculated from all sources in our final catalog, while the bottom 
panel shows w{0) calculated only from sources detected above a de- 
tection SNR threshold of 5. As expected, the results are consistent 
with zero correlation. 



7.5. Comparison with PAPER Results 

Comparing the prese nt results from th e MWA with re- 
sults from PAPER (Pa rsons et al.|[2010l ) is a particularly 
useful exercise, as both arrays are new, broadband, wide 
field-of-view instruments with similar uv coverage, and 
both are intended to be used to make EoR power spec- 
trum measurements. We find that 43 of the sources found 
in the survey (described in Section [5]) of iJacobs et al.l 
(|2011f) are located in our survey region. A search in our 
MWA-32T catalog reveals unique counterparts within 30' 



(the PAPER beam size used by IJacobs et all 120111 for 
source association) of the PAPER locations for 31 out of 
these 43 sources, multiple counterparts in 11 cases, and 
no counterpart in one case. 

The one source with no MWA counterpart is 24° from 
the center of the MWA Hydra A field, and corresponds to 
a local maximum in the MWA image; however it was not 
detected by the automatic source finding algorithm. For 
each of the 10 sources with multiple MWA counterparts, 
an estimate of the blended flux was obtained by sum- 
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Fig. 11. — Comparison of the fluxes presented in this work 
with those f r om th e 145 MHz PAPER source Hst presented in 
[Jacobs et al.l II2011I ). Error bars have been omitted from the plot 
for clarity. A total of 43 sources from the PAPER list are within 
the MWA field. MWA sources, which are within 30' of a PA- 
PER source are matched. A total of 31 of the PAPER sources 
had unique counterparts, while 11 PAPER sources matched with 
multiple MWA sources and 1 PAPER source did not have a de- 
tected MWA counterpart (although there is a local maximum in 
the MWA map at the location of the PAPER source). The MWA 
fluxes are calculated by summing the flux of all MWA sources which 
match with a PAPER source, and scaling the flux to a frequency 
of 145 MHz assuming a spectral index of a = —0.8 (5 <x !^~"). 
The solid fine shows the unity flux-ratio locus; the MWA sources 
are on average 17% brighter than the PAPER sources. This fltted 
flux ratio locus is plotted as a dashed Hne. 



ming the flux of all MWA sources within the PAPER 
beam. Other than Hydra A, all 41 PAPER fluxes are 
consistent with the MWA blended fluxes (the Hydra A 
flux reported in the PAPER catalog was corrupted by 
the filtering used in the PAPER analysis. [Jacobs! I2011L 
private communication). A weighted average of the ra- 
tios of the MWA and PAPER source fluxes yields the 
average ratio < S'mwa/S'paper > = 1.17 ±0.10. How- 
ever, we note that the PAPER flux scale was set using 
measurements of two calibration sources from the Culgo- 
ora source list, whereas the MWA-32T flux scale was set 
using a fit to an ense mble of Culgoora and MRC mea- 
surements. iSled (|1977D note that the Culgoora flux scale 
may be depressed by 10%, with additional flux uncer- 
tainties of between 13% and 39% for individual source 
measurements. If these Culgoora flux uncertainties are 
taken into account as potential errors on the PAPER flux 
scale, then the significance of the difference between the 
MWA and PAPER flux scales is decreased. 

The standard deviation of the MWA to PAPER flux ra- 
tios after correcting for the different flux scales is ~25%. 
This is smaller than our estimate of the MWA flux un- 
certainties based on flux predictions from the MRC and 
Culgoora measurements. This indicates that the flux 
comparison with the MRC and Culgoora lists may be 
affected to a considerable extent by radio source vari- 
ability or other systematic effects. 

7.6. Candidate Ultra Steep Spectrum Sources 



Fig. 12. — Distribution of spectral indices (5^ oc i/") be- 
tween 154.25 MHz and 4.85 GHz for sources identifled in this pa- 
per, based on a comparison with the Parkes-MIT-NRAO catalog 
l|Grifl5th &: Wright. 1993 ^. To avoid issues of source confusion and 
blending, only sources which could be unambiguously associated 
with single PMN counterparts are included in the histogram. 

"Ultra Steep Spectrum" (USS) radio sources form 
a comp elling class of candidat e high redshift radio 
sources jDe Brcuck ct'al. 2000; Di Matteo et al.' ^2003 
iDe Breuck^et al. 2002: Broderick et al. 2007 ) . Low fre- 
quency radio observa tions are par ticularly sensitive to 
these objects (see, e.g. lPedani|[2003() . We have conducted 
an analysis of the sources detected in the MWA flelds in 
an attempt to identify additional USS radio sources. A 
table of candidates is presented in Table [6l We calcu- 
lated spectral indices by using the PMN 4.85 GHz survey 
(jGriffith k Wright! [TOOl together with the MWA full- 
band average flux measurements. We associated MWA 
sources with PMN counterparts if their positions were co- 
incident within 5'. To avoid source confusion and blend- 
ing issues, we excluded any MWA sources with more than 
one PMN counterpart within a 30' radius. A total of 
331 sources were identified for which we could unambigu- 
ously identify a counterpart and extract a spectral index. 
A histogram of the spectral indices are shown in Fig- 
ure [121 This histogram appears consistent with the low- 
fre quency selected spectral index distributions obtained 
by (Pe Breuck et all (I2000D. and plotted in their Fig 



ure 7 (however, the IDe Breuck et all 120001 analysis used 
slightly different frequencies). Using a low- frequency se- 
lected distribution results in a sample which is signifi- 
cantly more sensitive to the USS sources. 

We choose a spectral index cutoff of a < —1.2 for USS 
source candidates, and find 3 sources which match the 
criteria. All three sources have counterparts in the MRC, 
and the source MWA J1032-3421 has a counterpart in 
both the PAPER and Culgoora source lists (although the 
other two USS candidates do not). A further 33 sources 
are identified which have no counterpart in the PMN cat- 
alog wit hin 1°. Using the PMN catalog limiting flux of 
50 mJy (jGriffith fc Wright I1993D for these sources, we 
find that 25 of these 33 sources have an inferred spectral 
index of a < —1.2. Of these 25 sources, 11 have unique 
counterparts in the MRC (one additional source has mul- 
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TABLE 6 

Ultra-Steep Spectrum Source Candidates 



Name 



RA 



DEC 



'MWA,avg 



"PMN 



■S'mrc 



"MRC 



Detection 
SNR Level 



J1009-1207 


10h09ml9s 


-12°07'46" 


10.69±3.21 0.17±0.01 -1.20±0.09 


3.23±0.11 


-1.23±0.31 


52.5 


J1032-3421 


10h32m60s 


-34° 21' 19" 


17.38±5.31 0.27±0.02 -1.21±0.09 


5.59±0.25 


-1.17±0.32 


15.5 


J1042-I-1201 


10h42m56s 


+12°01'30" 


15.55±4.80 


<-1.66±0.09 


8.90±0.37 


-0.57±0.32 


15.1 


J1034+1111 


10h34ml3s 


+11°11'23" 


5.89±1.92 


<-1.38±0.09 


3.80±0.16 


-0.45±0.34 


8.3 


JlOOO+1400 


10h00ml6s 


+14°00'17" 


7.54±2.51 


<-1.45±0.10 


3.06±0.10 


-0.93±0.34 


6.3 


J0831-2922 


08h31m24s 


-29°22'26" 


3.65±1.24 


<-1.24±0.10 


1.51±0.06 


-0.91±0.35 


5.7 


J1007+1246 


10h07m28s 


+12°46'50" 


7.46±2.52 


<-1.45±0.10 






5.6 


J0855+0552 


08h55ml8s 


+05°52'50" 


4.11±1.40 0.06±0.01 -1.23±0.12 


1.62±0.07 


-0.96±0.35 


5.4 


J0834-3443 


08h34m29s 


-34°43'22" 


6.14±2.12 


<-1.40±0.10 






5.0 


J0828-3201 


08h28ml4s 


-32"01'26" 


4.68±1.66 


<-1.32±0.10 






4.8 


JlOOl+1108 


lOhOlmOls 


+11°08'19" 


3.98±1.37 


<-1.27±0.10 


1.57±0.06 


-0.96±0.36 


4.6 


J0832-3326 


08h32m38s 


-33°26'00" 


5.14±1.84 


<-1.34±0.10 


2.30±0.09 


-0.83±0.37 


4.4 


J1008+1201 


10h08mlls 


+12°01'07" 


4.48±1.65 


<-1.30±0.11 






4.2 


J1028-I-1158 


10h28m32s 


+11°58'58" 


3.56±1.30 


<-1.24±0.11 






4.1 


J1112-I-1112 


Ilhl2m45s 


+11°12'45" 


4.88±1.92 


<-1.33±0.11 


1.45±0.07 


-1.25±0.41 


3.8 


J1021+1303 


10h21m34s 


+13°03'55" 


3.76±1.47 


<-1.25±0.11 






3.7 


J1034+1428 


10h34ml8s 


+14°28'48" 


4.68±1.85 


<-1.32±0.11 






3.5 


J1104+1103 


Ilh04m21s 


+11°03'31" 


5.35±2.01 


<-1.36±0.11 






3.4 


J0827-3322 


08h27m30s 


-33°22'43" 


3.61±1.47 


<-1.24±0.12 






3.4 


J0918+1226 


09hl8m48s 


+12°26'31" 


5.67±2.34 


<-1.37±0.12 


2.02±0.09 


-1.06±0.43 


3.3 


J1027+1347 


10h27ml0s 


+13°47'25" 


3.94±1.52 


<-1.27±0.11 






3.3 


J1114+1048 


Ilhl4m39s 


+10°48'59" 


4.28±1.71 


<-1.29±0.12 






3.2 


J0843+1115 


08h43m31s 


+11°15'20" 


6.34±2.43 


<-1.40±0.11 


1.34±0.07 


-1.60±0.40 


3.2 


J1015+1141 


10hl5m39s 


+11°41'05" 


3.16±1.31 


<-1.20±0.12 






3.2 


J0929+1133 


09h29m20s 


+11°33'39" 


5.71±2.13 


<-1.37±0.11 


1.69±0.08 


-1.25±0.39 


3.1 


J1019-I-1405 


10hl9m53s 


+14°05'41" 


3.69±1.46 


<-1.25±0.11 


0.87±0.06 


-1.49±0.41 


3.1 


J1026-I-1431 


10h26m53s 


+14°31'09" 


3.96±1.67 


<-1.27±0.12 






3.0 


J0918+1114 


09hl8m34s 


+11°14'45" 


5.40±2.24 


<-1.36±0.12 






3.0 



Note. — This USS sample is created by matching sources in the MWA-32T catalog with uniquely corresponding sources in the 
4.85 GHz PMN catalog. Sources which have spectral indices of a < —1.2 {S ex; i.*") are identified as ultra-steep candidates. The 
sources arc sorted in order of decreasing detection SNR level. MWA sources which are within 5' of a PMN source, and have no 
other PMN counterparts within 30 arc matched as counterparts and used to calculate the spectral index. Additionally, for MWA 
sources with no PMN source within 1°, a spectral index limit is calculated using the PMN flux limit of 50 mjy. For comparison, 
the ultra-steep candidates are matched with MRC candidates within 15' (sources with multiple matches within 30' are excluded), 
and an additional spectral index is calculated. 



tiple counterparts), 9 have counterparts in the Culgoora 
source list, and none have counterparts in the PAPER 
source hst. 

A comparison betw een the USS source hst of 
iDe Breuck et al.l ()2000f ) and the source hst identified 
in this work finds only 1 source from their list which 
matches a MWA USS candidate within 15': MWA J1032- 
3421. They select this source from an analysis of the 
MRC and PMN samples, and find a spectral index 
of -1.23 ± 0.04, which is consistent with the MWA- 
32T meas urements. There are an additional 21 sources 
from the IDe Breuck et alj (|2000| ) sample which match 
MWA sources, however only three of these MWA sources 
matched uniquely with a PMN counterpart: MWA 
Jl 133-2717, MWA J0941-1627, and MWA J0937-2243. 
These three source s were identified as US S sources in 
the "TN" sample of IDe Breuck et al.l (|2000l ). which used 
measurements at 365 MHz and 1.4 GHz, but they did 
not meet the USS criteria in the MWA-PMN compar- 
ison. The flatter MWA-PMN spectra may indicate a 
high-frequency turnover in th e source spectrum, similar 
to that noted for J0008-421 in lJacobs et all (|2011D . This 
interpretation is suppo rted by 74 MHz V LSS measure- 
ments of these sources (jCohen et al.ll200"7l ). which imply 
a spectral flattening at low frequencies, although further 
follow-up will be important to definitively establish this 



behavior. 

It is important to note that these MWA sources are 
only candidates, and should not be treated as definitive 
ultra-steep spectrum sources. Flux calibration and mea- 
surement errors as well as blending issues due to the low 
MWA resolution and time variability may result in a re- 
classification of these sources upon more detailed inves- 
tigation. Additionally any time variability or errors in- 
troduced by the different catalog resolutions and fitting 
algorithms will add to errors in this candidate list. The 
MWA instrument is under continued development, and 
the fidelity of these studies will improve as the system- 
atics are better understood. However, this candidate list 
serves as a good basis for more detailed follow-up and 
investigation. 

8. CONCLUSIONS AND FUTURE WORK 

The goals of this work were to verify the performance 
of the MWA subsystems and the MWA-32T system, to 
explore techniques for future EoR experiments, and to 
deepen our understanding of the radio sky at these fre- 
quencies. The analysis and results presented in this pa- 
per served to help commission MWA-32T and represent 
an assessment of its performance. Specifically, our abil- 
ity to successfully solve for antenna gains in spite of the 
wide field of view and direction-dependent primary beam 
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increases our confidence in our ability to calibrate the 
full 128-tile MWA array. The high level of agreement 
of our position measurements with existing source po- 
sition measurements provides further verification of our 
understanding of the geometry of the array and the cal- 
ibration procedures as well as our expectations for iono- 
spheric effects on scales relevant to the MWA-32T array. 
In short, the high fidelity, wide- field images produced in 
this analysis helps to build confidence that the MWA will 
be able to achieve its design goals. The measured fluxes 
in the maps agree with expectations from previous cat- 
alogs with a scatter of about 30%. The magnitude of 
these flux residuals is similar to what is reported from 
other pathfinder low-frequency arrays, however it is sig- 
nificantly larger than expected based on the noise in the 
maps. Further work is needed to understand the effects 
that cause these flux discrepancies. 

Future EoR experiments with the MWA will require 
long integrations, which will require the combination of 
data from many pointings. We have shown in this work 
that these data can be corrected for primary beam ef- 
fects and weighted averages can be formed to increase the 
sensitivity and fidelity of the resulting maps. In this ini- 
tial exploration we made certain simplifying assumptions 
that allowed us to do this analysis in the image domain. 
Future, deeper, investigations will require more sophis- 
ticated techniques that account for differences between 
individual antenna elements. Developing and verifying 
these techniques should also be a priority. 

A distinguishing feature of the images presented in this 
paper is the degree to which they are confusion-limited. 
The effects of source confusion are not expected to be an 
issue for EoR experiments. Simulations indicate, though 
it remains to be shown in practice, that subtracting 
the brightest sources, and treating the fainter sources as 
smooth contributors to each pixel will allow an adequate 
separation of foregrou nds from the EoR signal (see, e.g. 
iLiu fc Tegmark|[2011[ ). For other science goals, however, 
it is quite clear that for the low-frequency arrays planned 
for the near future, the confusion limit of continuum im- 
ages will be reached in relatively short integration times. 
Thus, much of the science will have to contend with the 
difficulties of measurements in crowded fields. We have 
presented the results of an automatic source extraction 
algorithm. While the results are quite good, with relia- 
bility over 90%, further development of algorithms that 
extract more accurate source model s from crowded fields 
(such as the algorithm presented in lHancock et al.ll2012| ) 
should be a priority. 

The source list presented in this paper serves as the 
deepest catalog of radio sources generated from a blind 
search in this region of the southern sky in this frequency 
band and is the first MWA characterization of the EoR 
point source foreground in the MWA EoR field. Further 
refinement of statistical descriptions of this component 
of the foreground should be possible, and represent a 
promising avenue for future work. The production of this 



large, high quality survey with only ~25 hours of data 
highlights the survey power of this instrument. The tech- 
niques used to survey these two fields can be extended to 
complete an MWA all-sky survey. It is important to note 
that because this survey was carried out near the source 
confusion limit, many of the fainter sources in the sam- 
ple are likely blends of multiple sources. The USS can- 
didates identified in this work, although likely affected 
by this blending, serve as a set of candidates for high- 
redshift radio sources, and are good targets for follow-up 
with higher resolution low-frequency instruments. 

The MWA is in the process of a buildout to a 128- 
tile (128T) interferometer. The 128T array will have 
four times the number of collecting elements as the 32T 
array, and the maximum baseline length will be increased 
to ~3 km. The higher resolution of the 128T array will 
yield maps with a lower source confusion limit, and will 
enable a deeper survey of this sky field with a reduced 
number of blended sources. With the added sensitivity, 
the 128T array will still be able to rapidly produce source 
confusion limited maps. The techniques developed as 
part of the work described herein will allow evaluation 
of the 128T instrument as it is commissioned. As the 
quality of the sky model at these frequencies improves, 
the increased calibration accuracy and ability to subtract 
foregrounds from the data will move us closer towards the 
goal of detecting the 21-cm signal during reionization. 
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